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j^ . Abstract 

t^ ■ The Pauli exclusion principle is advocated for constructing the proton and 

neutron deep inelastic structure functions in terms of Fermi-Dirac distribu- 
O ■ tions that we parametrize with very few parameters. It allows a fair descrip- 

^ ; tion of the recent NMC data on Fi{x, Q^) and F:^{x, Q^) at Q^ = 4,GeV'^, as 

well as the CCFR neutrino data at Q^ = 3 and 5GeV^. We also make some 



p i' reasonable and simple assumptions to relate unpolarized and polarized quark 



parton distributions and we obtain, with no additional free parameters, the 



(D . spin dependent structure functions xg^{x, Q^) and xg^{x, Q^). Using the cor- 

' rect Q^ evolution, we have checked that they are in excellent agreement with 

the very recent SMC proton data at Q^ = lOGeV^ and the SLAC neutron 
X ■ data at Q^ = 2GeV\ 
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Many years ago Feynman and Field made the conjecture^^^ that the quark 
sea in the proton may not be flavor symmetric, more precisely li > m, as a 
consequence of Pauli principle which favors dd pairs with respect to uu pairs 
because of the presence of two valence u quarks and only one valence d 
quark in the proton. This idea was confirmed by the results of the NMC 
experiment [^1 on the measurement of proton and neutron unpolarized struc- 
ture functions, -^2(0;). It yields a fair evidence for a defect in the Gottfried 
sum rule'^1 and one finds 

I^= I — [F|(x) - F^{x)] = 0.235 ± 0.026 (1) 

instead of the value 1/3 predicted with a flavor symmetric sea, since we have 
in fact 

lG = -{u + u-d-d) = - + -iu-d). (2) 

A crucial role of Pauli principle may also be advocated to explain the well 
known dominance of u over d quarks at high a;,^ which explains the rapid 
decrease of the ratio F2{x)/F2{x) in this region. Let us denote by q^{q^), u 
or d quarks with helicity parallel (antiparallel) to the proton helicity. The 
double helicity asymmetry measured in polarized muon (electron) - polarized 
proton deep inelastic scattering allows the determination of the quantity 
^{x) which increases towards one for high x,'^'^' suggesting that in this 
region u' dominates over m^, a fortiori dominates over d^ and d^, and we will 
see now, how it is possible to make these considerations more quantitative. 
Indeed at Q^ = the first moments of the valence quarks are related to the 
values of the axial couplings 



t I ,t l + F -D ,1 I- F + D , , 
«L = 1 + F, «ii = l-F, dl, = ^ , dl, = ^ , (3) 

so by taking F = 1/2 and D = 3/4 (rather near to the quoted valuesf^' 
0.461 ± 0.014 and 0.798 ± 0.013) one has ul^^ = 3/2 and m^^j = 1/2 which 
is at the center of the rather narrow range (c^vai' '^vai) ~ (3/8,5/8). The 
abundance of each of these four valence quark species, denoted by Pvai, is 
given by eq. (pi) and we assume that the distributions at high Q^ "keep a 



memory" of the properties of the valence quarks, which is reasonable since for 
X > 0.2 the sea is rather small. So we may write for the parton distributions 

p(x) = F(x,pvai) (4) 

where F is an increasing function of Pvai- The fact that the dominant distri- 
bution at high X is just the one corresponding to the highest value of p^ai, 
gives the correlation abundance - shape suggested by Pauli principle, so we 
expect broader shapes for more abundant partons. If F(a;,Pvai) is a smooth 
function of Pvai, its value at the center of a narrow range is given, to a good 
approximation, by half the sum of the values at the extrema, which then 
implies^^l 



This leads to 



4ai(a^) = 2^vai(a;). (5) 



Au^^i{x) = ul^i{x) - ui^i{x) = Mvai(a;) - d^s^iix) (6) 



and, in order to generalize this relation to the whole u quark distribution, we 
assume that eq. (|^) should also hold for quark sea and antiquark distributions, 
so we have 

^Usea{x) = Au{x) = u{x) — d{x) . (7) 

Moreover as a natural consequence of eq. (§), we will assumeQ 

Ad,,i{x) = (F - D)dUx) . (8) 

Finally we will suppose that the d sea quarks (and antiquarks) are not po- 
larized i.e. 

A4ea(a:) = Ad{x) = (9) 

and similarly for the strange quarks 

As{x) = As{x) = . (10) 



^It is amusing to remark that with the values of F and D quoted above, we have in 
fact Advai{x) = —hdvaiix) which coincides with the so-called conservative SU{6) model^^l! 



Clearly the above simple relations (6)- (10) are enough for fixing the deter- 
mination of the spin dependent structure functions xgi'"'{x,Q'^), in terms 
of the spin average quark parton distributions. We now proceed to present 
our approach for constructing the nucleon structure functions ^^'"(x, Q^), 
xF^'^{x, Q"^), etc... in terms of Fermi-Dirac distributions which is motivated 
by the importance of the Pauli exclusion principle, as we stressed above. A 
first attempt for such a construction was made in ref. [10], but here, as we 
shall see, our method is slightly different and leads to significant improve- 
ments. Let us consider u quarks and antiquarks only, and let us assume 
that at fixed Q'^,ul^i{x), ul^i{x), u^{x) and U'^{x) are expressed in terms of 
Fermi-Dirac distributions, in the scaling variable x, of the form 

xp(x) = apX^^ / {exp{{x — x{p))/x) + 1) . (11) 

Here x{p) plays the role of the "thermodynamical potential" for the fermionic 
parton p and x is the "temperature" which is a universal constant. Since va- 
lence quarks and sea quarks have very different x dependences, we expect 
< 6p < 1 for u\'.^x{x) and 6p < for v}'^{x). Moreover x{p) is a con- 
stant for Mvai(^)' whereas for u^'^(x), it has a smooth x dependence. This 
might reflects, the fact that parton distributions contain two phases, a gas 
contributing to the non singlet part with a constant potential and a liquid, 
which prevails at low x, contributing to the singlet part with a potential 
slowly varying in x, that we take linear in ^/x. In addition, in a statistical 
model of the nucleonl^^l , we expect quarks and antiquarks to have opposite 
potentials, consequently the gluon, which produces qq pairs, will have a zero 
potential. Moreover since in the process G — > (jsea + ^, Q'sea and q have oppo- 
site heli cities, we expect the potentials for ul^^ (or v}) and u^ (or u^^^ to be 
opposite. So we take 

x{v)) = — x(n^) = xq + X\\fx . (12) 

The d quarks and antiquarks are obtained by using eqs. (|) and (|^) and 
concerning the strange quarks, we take in accordance with the data '^^1 

u[x) + d[x) 
s[x) = s{x) = . (13) 

Finally for the gluon distribution, for the sake of consistency, one should 



assume a Bose-Einstein expression given by 

-G(^) = ^f^, (14) 

with the same temperature x and a vanishing potential, as we discussed 
above. 

To determine our parameters we have used the most recent NMC data^^^ 
on F2{x) and -F^(x) at Q^ = AGeV"^ together with the most accurate neutrino 
data from CCFRt^^'^^J on xF^^(x) and the antiquark distribution xq{x) t^^'. 

The universal temperature is found to be 

X = 0.120 (15) 

and for valence quarks we get the three free parameters 

b{ul,) = l/2b{uij = 0.399, x{ul,) = 0.502, i{ul,) = 0.163. (16) 

This relation between the b's is imposed by the small x behavior of xF^^{x), 
a^ and a^ are not free parameters, but two normalization constants which are 
fixed from the obvious requierements to have the correct number of valence 
quarks in the proton. As we noticed before ul^i{x) dominates, so it is not 
surprising to find that it has a larger potential than ul^■^{x) 0. 
For antiquarks we have four additional free parameters 

b = -0.358, a^ = 0.024, xq = 0.215 and xi = -0.388 for ul (17) 

b is the same for u^ and u^. 

When a; — >■ 0, from Pomeron universality, one expects xu{x) = xd{x) ^ 0, 
so d^ is determined by this constraint. 

We show the results of our fit for -Ff (x) — -F^f (x) and ^^(a;)/^^^) by the 
solid lines in Figs, la and b and for xF^'^lx) and xq{x) in Figs. 2a and b. 
The accuracy of these neutrino data gives strong constraints on both valence 
and sea quark distributions. The description of the data is very satisfactory 
taking into account the fact that we only have eight free parameters and this 
certainly speaks for Fermi-Dirac distributions. Note that we find Iq = 0.232 



■^In a statistical model of the nucleoli!^ ^1, the potentials associated with u and d quarks 
are taken in the ratio 2-^^^ which is a much smaller than the value of ai(^tvai)/^("vai) ~ "^ 
we have found. 



in beautiful agreement with eq. ([1|). The steady rise oi xq{x) at small x leads 
to a rise of Ff which is consistent with the first results from Hera'^^'. For the 
fraction of the total momentum carried by quarks and antiquarks we find 

/o xu{x)dx = 0.298, /q xd{x)dx = 0.115, (18) 

Jq x[u{x) + d{x) + s{x) + s{x)]dx = 0.147. 

Concerning the gluon distribution, we find ac = 14.536 and bo = 0.912 
and xG{x) is fairly consistent with some preliminary indirect experimental 
determination from direct photon productiont^^^ , from neutrino deep inelastic 
scattering t^^l at Q^ = 5GeV^ and at high Q^ and smaller x from NMC'"^''''. 

Let us now turn to the polarized structure functions xgi'"'{x, Q^) which 
will allow to test our simple relations (|6D-(|T0|). Since the proton data from 
EMCl^l and SMCt^l are at Q^ = lOGeV^ and the neutron data from SLACl^^lQ 
is at Q^ = 2GeV^, we have to consider the Q^ evolution in order to use our 
parton distributions determined at Q^ = AGeV'^. For this purpose we have 
used a numerical solutionl^^l of the Altarelli-Parisi equations^^^' which lead to 
relatively small corrections in the Q^ range we are dealing with. In Figs, la, b 
the dashed lines are the theoretical predictions at Q^ = lOGeV^. As expected 
we see that for low x, F2{x) — -F^(x) increases with Q^ whereas it decreases 
with Q^ for high x, leaving the integral unchanged. The Q^ dependence of 
the ratio -F^(x)/F2'(x) has the right trend although probably a bit too weak 
compared to experimental observation which has been attributed to different 
higher twist effects for proton and neutron'^^l 

At this stage we would like to examine the consequence, of our simple re- 
lations (|)-(|iy). If the d (valence and sea) quarks were unpolarized, eqs. (|^) 
and (0) allow to relate the contribution of u quarks to xgi{x), to the contri- 
butions of u and d to -^^(a;) — -F;f (x) i.e. 

X9'i{x)\u = I (F^ix) - F^{x)) U, . (19) 

First, by comparing Fig. la and Fig. 3 Q one sees very clearly, the similarity 
of the two sets of data pointsQ. 

•^There is also some SMC data'^^^ on the polarized structure function xgf on deuterium. 

''The vertical scales have been chosen in such a way that one absorbs the factor 2/3 by 
superimposing the two figures. 

^This was first noticed in ref. [8] but, with more accurate data, it becomes now very 
convincing. 



Second, on Fig. 3 the dashed hne represents simply the dashed hne of 
Fig. la multiphed by a factor 2/3, whereas the dotted hne corresponds to 
the case of a flavor symmetric sea, i.e. u{x) = d{x), or, in other words, to a 
zero polarization of the u sea quarks. This shows why we strongly suspect 
that the defect in the Gottfried sum rule and the defect in the proton EUis- 
Jaffe sum rule'^^' are closely related. We still think it has nothing to do with 
the polarization of the strange quarks, that we took to be zero (see eq. ([10|)), 
which is supported by reasonable phenomenological argumentsf^^^ . Moreover, 
in this approach, the strange quarks do not even participate, because they 
cancel in the difference F2{x) — -F^(x). Finally, if one takes into account the 
polarization of the d valence quarks by using eq. (||), we get the solid hne 
in Fig. 3 which improves the agreement with the data. In fact we found for 

rO.J 

/ g^^{x)dx = 0.136 (20) 

J0.003 

in beautiful agreement with eq. (^) of ref. [6]. 

Concerning the neutron polarized structure function xg^lx) we show in 
Fig. 4 a comparison of the SLAC data'^^^ at Q^ = 2GeV^ with our theoretical 
calculations. The dashed hne corresponds to the case where d quarks are 
assumed to be unpolarized and it clearly disagrees with the data. However 
by including the d valence quark polarization according to eq. (P), we obtain 
the solid line in perfect agreement with the data and we find for Q^ = 2GeV^ 

1 
g'^{x)dx = -0.017 . (21) 



To summarize we have given an accurate description of deep inelastic 
scattering data at low Q^ in terms of Fermi-Dirac distributions parametrized 
with only eight free parameters for quarks and antiquarks. Although we 
have some understanding of their meaning, much remains to be done for a 
more fundamental theoretical interpretation, in terms of new information for 
the nucleon structure. We have proposed a set of simple relations between 
unpolarized and polarized quark (antiquark) distributions for which, so far, 
there is a striking experimental evidence. Of course our approach has to 
be further tested with more accurate deep inelastic scattering data and in 
particular the important issue of the validity of the Bjorken sum rulef^^^. 
Polarized proton collisions at high energies will also provide independent 
tests which will be most welcome in the near future'^^'. 
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Figure Captions 

Fig. la The difference F2{x) — -F^(a;) at Q^ = AGeV"^ versus x. Data are from 
ref.[2] and the solid hne is the result of our fit. The dashed hne is the 
theoretical result after evolution at Q^ = lOGeV^. 

Fig. lb The ratio F2"(a;)/F|'(a;) at Q^ = AGeV"^ versus x. Data are from ref.[2] 
and the solid line is the result of our fit. The dashed line is the theo- 
retical result after evolution at Q^ = lOGeV^. 

Fig.2a The structure function xF^'^{x) versus x. Data are from ref. [13] at 
Q^ = 3GeV^ and the solid line is the result of our fit. 

Fig.2b The antiquark contribution xq{x) = xu{x) + xd{x) + xs{x) at Q^ = 
3GeV^ (full circles) and Q^ = bGeV"^ (full triangles) versus x. Data 
are from ref. [12] and solid line is the result of our fit. 

Fig. 3 xgi{x) at < Q^ >= lOGeV^ versus x. Data are from ref. [5] (full 
squares) and ref. [6] (full circles) together with our predictions at Q^ = 
lOGeV^. (Dotted line is the contribution of Au^gx{x) only, dashed line is 
the contribution of Au{x) and Au{x), and solid line is the contribution 
of Au{x), Au{x) and Ad^3x{x)). 

Fig. 4 xg'^{x) at < Q^ >= 2GeV^ versus x. Data are from ref. [18] together 
with our predictions at Q^ = 2GeV^ (Dashed line is the contribution 
of Au{x) and Au{x) only and solid line is the contribution of Au{x), 
Au{x) and Ad^aiix)). 
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